Dynamical local field, compressibility and frequency sum rules for quasiparticles 
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The finite temperature dynamical response function including the dynamical local field is derived 
within a quasiparticle picture for interacting one-, two- and three dimensional Fermi systems. The 
correlations are assumed to be given by a density dependent effective mass, quasiparticle energy shift 
and relaxation time. The latter one describes disorder or coUisional effects. This parameterization 
of correlations includes local density functionals as a special case and is therefore applicable for 
density functional theories. With a single static local field, the third order frequency sum rule 
can be fulfilled simultaneously with the compressibility sum rule by relating the effective mass and 
quasiparticle energy shift to the structure function or pair correlation function. Consequently, solely 
local density functionals without taking into account effective masses cannot fulfill both sum rules 
simultaneously with a static local field. The comparison to the Monte-Carlo data seems to support 
such quasiparticle picture. 
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I. INTRODUCTION 

The response of an interacting Fermi system with the 
potential Vq to an external perturbation is the basic 
source of our knowledge about the interaction and dy- 
namical as well as statical properties of the system. This 
response function has been therefore a central issue of 
many body theories. 

The density response function gives the variation of 
the density in terms of the external potential 



(1) 



The polarization is defined as the density variation in 
terms of the induced potential 



Sn{q,Lo) = niq,Lu)SV"'''{q,uj) 



(2) 



where we suppress the notation of obvious g-dependence 
in the following. The induced potential itself is the sum 
of the external potential and the effective interaction po- 
tential {Vq -\- fq{uj))Sn 



6V"'%Lo) = [Vq + fq{uj))5n{uj) + V^'iio). 
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Therefore, from (||)-(||) we have the relation between re- 
sponse and polarization 



l-{Vq + fq{0j))Ii{iOy 



(4) 



The local field fq{Lo) describes the modification in the 
restoring force brought about by particle correlations. 
This field prevents the particles frorrusampling the full 
effect of interaction at short distanceail. 

The dielectric function relates now the induced densi- 
ties to the external potential via 
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Vq6n{uj) 
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such that the dielectric function reads 



Vqn{u) 



1 - fq{uj)n{L0) 
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The theoretical effort consists in determining the lo- 
cal field fq{Lo) which represents the local correlation 
and which depleted the induced potential by fq{uj) = 
—G{u!)Vq. As long as this local field is a dynamical one 
this is an exact relation. The different theoretical treat- 
ments differ in this local field corrections, for an overview 
seecl. Mostly static approximations, fq = fq{0) — —VqG, 
have been proposed in the nast. ft has started with the 
pioneering work of Hubbardtl who first introduced the no- 
tation of local field and took into account the exchange- 
hole correction resulting in 



G 



H 



1 g2 



1 9' 
2fe| 



■ + oiq"^) 
o{l/q^) 



(7) 



While this expression has established a remarkable im- 
provement of the dielectric function in random phase ap- 
proximation (RPA), it has been soon recognized insufh- 
cient due to the lack of self-consistency which leads the 
pair correlation function still to unphysical neaative val- 
ues. This has been repaired by Singwi et. al.u by using 
exchange-correlations 
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with 



(9) 



f 



where s is the spin degeneracy and where the static struc- 
ture factor 
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with inverse temperature f3 = l/T is hnked to the pair 
correlation function via 



ffr - 1 = 



1 



dq 



■e'""^(5, - 1). 
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This provides a self-consistent problem in solving dielec- 
tric function, structure function and static local field si- 
multaneously. The advantages of this result compared to 
the Hubbard result with respect to the pair correlation 
function and large wavevector limit has been discussed 
bytl. Recent comparisons with molecular dynamics sim- 
ulations for a hard sphere gas is presented inQ where a 
good agreement is found for thermodynamical properties. 

The expression (||) Jias been improved further by 
Pathak and VashishtaQ demanding that the responsg 
function should fulfill the third order frequency sum ruleEl 
which resulted into 
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leading to the improved small distance limit discussed 
irJa. The difference at short distance to Gstls comes 
from the motion of particles inside the correkttion hole 
which is condensed in the dynamical behaviouij. 

If one takes into account the differeiiGe-|between un- 
correlated and correlated kinetic energyliirS one obtains 
an additional —/S.E = —:^^y-{Eint — E) term to Gpv- 
This comes from the difference in correlated and uncor- 
related occupation numbers which can be expressed by a 
coupling constant integration and can be linked via the 
virial thcorcna_jta density derivatives of the pair correla- 
tion functioi£3lij. 

Parallel to the above discussions there has been dif- 
ferent iMpcpuements to derive local fields from the virial 
formulaa^tlHU which have resulted into expressions known 
from density variations 
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d"^ Eq / d{e'^)'^ < 0. The same conclusion are obtained inil 
using the virial theorem. 

Therefore the concentration is now mostly fo- 
cussed onj-tbe j-construction of dynamical local field 
correctional3't3ll3. The quantum versions of the Singwi- 
Tosi-Land-Sjolander, Gstls, aa d Vah ishta-Singwi, Gvs, 
theories have been discussed intSHljtd. These lead to pos- 
itive values for the pair distribution jfunction at short dis- 
tances valid for rather low densitiest£l. There it has been 
focussed on dynamical properties of the dynamic local 
field. While the high frequency limit is monotonic and 
similar to STLS and VS, the static limit can even exhibit 
peaked structures which can give rise to charge density 
waves underlying the nontrivial character of dynamical 
behavior. Unfortunately even the dynamic quantum ver- 
sion of the Singwi-Tosi-Land-Sjolapder local field cannot 
fulfill the compressibility sum ruleEl. We will show here 
that from a dynamical local field one can derive a static 
local field fulfilling both sum rules simultaneously if one 
takes into account the effective mass. This will resolve 
the puzzle of sum rules. 

Recent improvements of the response functionr-Skre-ba- 
sically due to numerical studies [af|-Monte CarlcO" £3 or 
molecular dynamical simulationsE£ll3. An interesting first 
principle numerical scheme is to solve the time depen- 
dent^adanoff and Baym equations including an external 
fieldED. Due to the variation of internal lines, already a 
Born diagram leads to a linear response which includes 
high order vortex corrections fulfilling sum rules consis- 
tently. The third order frequency sum rule regains imnpr r 
tance for reduced dimensional layered structurea2a'c3LJ. 
All results in this paper here can be straight forwardly 
generalized also to reduced dimensions as given in ap- 
pendix ^ for one, two and three dimensions. This could 
have an impacjtjjn recent discussions of two-layered elec- 
tron gassesE2ic3. 

Here we want to return to the analytical investiga- 
tions and will show that there exist a possibility to fulfill 
with one static local field correction both requirements, 
the third order sum rule and the compressibility sum 
rule. This is performed by working within a quasiparti- 
cle picture determining the effective mass appropriately. 
Within the frame of the quasiparticle picture we will de- 
rive an explicit expression for the dynamical local field 
factor which leads to the desired static limits. We obtain 
the identity 
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with various 1/2 < a < 1, seeEj. This procedure satisfy 
the eopipressibility sum rule almost exactly. 

InEl it was shown for Coulomb systems that one can- 
not construct a static local field factor which fulfills 
both limits, the compressibility and the third order sum 
rule (B2) since it would violate the theorem of Ferrell, 



with moments of the polarization given in appendix ^ 
and the effective mass m(n) and selfenergy A(n). It will 
provide a recipe how to construct a quasiparticle pic- 
ture by the knowledge of the structure factaiL.a.t small 
distances from experiments or simulationsrnrTH. This 
in turn leads to an easy microscopic parameterization in 
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terms of the effective mass and quasiparticle energy shift 
which could be compared directly to microscopic theo- 
ries. 

The unde rlyjnj g principle is analogous to the one found 
in literatureBEJ where a response function was param- 
eterized explicitly fulfilling sum rules and different con- 
straints. Other parameterizations can be found inE3 from 
variational approaches which are exact in the high den- 
sity limit. A different line of constructing the response 
function uses tiin^equency moments resulting in recur- 
rence relationaSa-Ea. 

Here in this paper we will give an alternative approach 
which uses general parameterizations of the selfenergy in 
terms of a functional which might depend on the density, 
energy and current. We restrict here to a one component 
system though the generalization to multicomponent sys- 
tems is straight, forwardLZIES and considered in different 
approacheg^jcll. 

In the next chapter we review shortly the compress- 
ibility sum rule and the third order frequency sum rule. 
In chapter III we give the dynamical response for quasi- 
particles which-is a special case of the general structure 
derived earliercSl. We show that the correct compressibil- 
ity appears and the third order sum rule can be satisfied 
if the effective mass and quasiparticle energy is chosen 
appropriately. Chapter IV will present some numerical 
results on the unpolarized electron gas at zero tempera- 
ture and the comparison with Monte Carlo simulations 
are discussed. Chapter V summarizes the results and in 
appendix^ we give frequently occurring correlation func- 
tions and practical forms for calculation in one-, two- 
and three dimensions. While all formulas in the main 
text are written as three dimensional, they hold for one- 
and two dimensions as well. Only the explicit correla- 
tion functions in appendix ^ have to be used as outlined 
there. Also the often required long wavelength expan sion 
of these correlation functions are given in appendix Bl . 
Appendix |c| finally is devoted to the short sketch of per- 
turbation theory and the derivation of the used sum rules 
for one, two and three dimensions. 



II. DETERMINATION OF STATIC LOCAL FIELD 
FACTOR 

Let us discuss two different boundaries for the static 
limit of the local field. This will be the compressibility 
and the third order frequency sum rule. 



A. Compressibility 

First we have to know how the compressibility should 
look like. This is particularly simple in the quasiparticle 
picture which we will use. In the quasiparticle picture the 
one-particle distribution function is a Fermi distribution 



where the density dependent effective mass m and the 
selfenergy shift A are obtained either from microscopic 
calculations or, as proposed here, from the sum rules. 
Thorough the paper we will understand now the masses 
as effective masses. 

From microscopic approaches the effective selfenergy 
shift and the effective mass are coming from the knowl- 
edge of the selfenergy a(ji, to) which determines the quasi- 
particle energy e via 
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2mQ 



cr(p, e). 



(16) 



The velocity of the quasiparticles is given by de/dp which 
leads to the definition of the effective mass 
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where the momenta and energy are set to the Fermi mo- 
menta and energy after derivatives. Consequently the 

2 

quasiparticle energy can be approximated by e « + A 
with the effective mass ([l7| ) and the energy shift A — 
a{pf,ef). This approximation has to be replaced by ther- 
mal averaging when finite temperature systems are con- 
sidered. A useful method would also be to read off the 
quasiparticle parameterization from c yirjei it parameteri- 
zations of the momentum distributionHcj. 

Using the definition of the compressibility one obtains 
directly from (R^ 
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where the free compressibility is = ^ j j^^F{p){l — 
F{p)). Alternatively, the energy shift A and the effective 
mass can be expressed by Landau parameter. 

In the next chapter we will present a consistent dynam- 
ical response function for the quasiparticle picture such 
that the correlations are parameterized by quasiparticles 
with an effective mass, an energy shift and a relaxation 
time. From this we will obtain the correct compressibil- 
ity ( p^ from the dynamical response via the static limit 
obeying the frequency sum rules. Actually, a static local 
field can be constructed provided we choose the effective 
mass appropriately. This will lead to a recipe how the 
effective mass can be determined from the structure fac- 
tor which is well known from Monte-Carlo simulations or 
experiments. _ 

The conventional compressibility sum ruleo reads with 
lim n = -n^/Co and ^ 



F{p) - (e 
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lim / Inie(w') = lim Re „ f^N-^ll r 



lim 

5-0 1 + /,(0)n2/Co 

(19) 



such that we can expect from the correct result ( p^ ) that 
the static local field has the form 



lim/,(0) = — 
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We will present a dynamical local field which leads in the 
static limit exactly to this desired result (pT 



where I{q) is usually presented in literaturi 
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and Sk is the structure factor (|lO|). 

In order to unders tand the different contributions, the 
short distance limit (C26) from ( ^6|) is performed 
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B. Frequency sum rules 

The sum rules can be easily read off from the fact that 
the response function is an analytical function in the up- 
per half plane and falls off with large frequencies faster 
than l/cj^ such that the compact Kramers Kronig rela- 
tion reads 



uj' — uj -\- i\j 



(21) 



closing the contour of integration in the upper half plane. 
From this one has 



Rex(cj) = 



with the moments 



duj' Imx(ij-'') 
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duj 



(22) 



(23) 



The first two moments are known exactly to be (ap- 
pendix ^ 



dtu nq 
< Lo >= I — ujiiRxy^) — 
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with-i-jtie density n and the mass to of the particles 
andEl&l 



J 



4 6 

q nq° 



^2 4j^3 



9 4- 

n'^q^ 



Vqiiq). 



(25) 



Here Emt is the kinetic energy of the interacting system 
and 

w ^ I f dk Ak-q?Vk 
I{q)=-- I j^^^{bk-q-i>k + ndk^q-ndkfi)- 



n J (27r)3 
I{q) - 1 



q^ Vq 
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where the last term comes from the 5^^ terpjuJIhe first 
term alone is sometimes called exact resultd'Oalld, which 
holds only for static local fields. The second term de- 
scribes the motion of particles inside the correlation hole 
and takes for X^oulomb just ~ go) which has been 

pointed out ira. Together one obtains the small distance 
result 



lim /(g)= (l-go)-l 



(29) 



in agreement with the direct expansion ( |C2q ). 

Now we proceed and derive the boundaries for the local 
field fq{uj) from (^ by the above sum rules. Therefore 
we look at the large uj expansion of (^) from which we 
can check with the help of ( p2|) the desired sum rulcsj 2^) 
and (|2|). The simple RPA leads to [see also (|^-(^] 
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+ o{l/u'>), (30) 



from which one gets with 



x(^) 



nq 



2E- 



+o(l/w'). 



nq^ n^q\,, , , ,,\ 1 
^ + ^(^g + /g(00))j-, 

(31) 



The first order energy weighted sum rule (24) is fulfilled 



trivially provided 11 fulfills it. The third order sum rule 
( p5| ) (Wi be fulfilled if we construct the local field accord- 
ing tc 



/,(oo) = -V,{1 + /(g)) - —{E - E„,t) 
= -VqI{q)-\{E-E,nt) 



(32) 



where E is the kinetic energy of the noninteracting sys- 
tem. The last term describes the fact that the third order 
frequency sum rule of the polarization function yields the 
noninteracting kinetic energy. iThis form neglecting the 
last term has been discussed inQ. In the later derivation 



4 



of the polarization function we cannot consider the ki- 
netic energy anymore as interaction free, since the relax- 
ation time appears as well as the effective mass. There- 
fore within the quasiparticle picture used here the differ- 
ence E — Eint vanishes or positively stated, is accounted 
for by the effective mass. To facilitate the comparison 
with the literature we have kept this difference formal as 

Ai? = —yj{E — Eint)- 

With (p2) we have given the constraint on the dynam- 
ical local field from the third order frequency sum rule. 
In the following we will present a dynamical local field 
which fulfills both requirements, the compressibility ( pO| ) 
and the frequency sum rule (p2h. 



III. DYNAMICAL RESPONSE FUNCTION 

was given the polarization function for an inter- 
acting quantum system imposing conservation laws on 
the relaxation time approximation. These polarization 
functions we have denoted by 11" for density conserva- 
tion imposed, II"'-' for density and current conservation 
and n"'-''^ for density, current and energy conservation. 
In the former paper we could give only formal matrix 
expressions for the response function. In appendix ^ we 
repeat shortly the way of derivation from the quantum 
kinetic theory and give now the e xplic it for m of the re- 
sponse function. We obtain with ( A15 ) and ( A16 ) 



1 - VoW'i-^{uj) - 2mt4ni3(w) 
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where 
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and 



Vo = Ia-v4 
on 4 
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The explicit expressions of the moments of the polariza- 
tion function are summarized in appendix ^ for compu- 
tation. From the response function ( ^3| ) we will read off 
the main result of this paper: the dynamical local field for 
quasiparticles with the effective mass and energy. Now 
we are going to work out explicitly the form of local field 
to show that the third order sum rule can be fulfilled and 
the correct compressibility is obtained. 



A. Dynamical local field 



Comparing (p3|) with (m an intermediate dynamical 



local field can be read off as [Vq 
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where we used the expansio n of appen dix H , explicitly 
( B16 ), in the first line and ( lB14| ) and ( |B15^ in the last 

1 ii-u-i 1 o r-1 v/^m i-\ -vlv- 4-1-1 ri 4- I I ^1 1 .1 i -t-J i + Oi^l-f- <^ vi + t-i m-i o /-iz-m^t-z-iI n 4- 1 c 



line. Please remark that 11"'^' itself contains correlations 
beyond the polarization in RPA, IIq. This we will present 
in a moment. 

First, we see the astonishing result that obviously 



lim fq{uj) = lim fg{uj) 



(37) 



the static local field required for the compressibility 
agrees with the infinite frequency limit required for the 
third order sum rule (^ij) . ^ This shows that the answer 
to the sumrule puzzle is not due to the explicit dynam- 
ical character of the local field as often claimed in the 
literature. Instead we will see in paragraph HID that it 



is due to the underlying selfenergy correction which has 
to obey certain relations. 

At this point it is important to avoid a misunderstand- 
ing. The intermediate dynamical local field, fg, is not 
the total one describing correlations beyond RPA polar- 
ization function, Xlg, which would be the case only in the 
infinite frequency limit. Instead, part of the correlations 
are already captured in Il'^'j'^. To make this explicit we 
write (BS 



1 



1 



(38) 



where the difference between 11"'^'^ aiid, IIq has been 
recasted into a local field contributions derived from 



(|A12|)-(|A14D 



*One should not be mislead to the conclusion that this vio- 
lates the Kramers Kronig relation for fq{uj). A toy example 
of Im/ — sin (aaj)/(tj^ — 47r^/a^) and the corresponding real 
part obtained from the Kramers Kronig relation (|2^ ) shows 
that indeed (^) can hold simultaneously with the Kramers 
Kronig relation. 
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with the last hne vahd for zero temperature. Oppositely 
in the static hmit /*(0) = 0. Together with (H) we 
obtain an effective local field renormalizing the RPA 
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where the last line is again the zero temperature limit. 
The high frequency limit required for the third order sum 
rule all agree /°*(oo) = fq{oo) = fq{oo). 



B. Connection to density functionals 

In order to establish the connection to the ground state 
exact relationsLJ of exchange correlation energy, £xc) we 
see from (40) 
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with ff{Q) < /^"(oo) < 0. This establishes-|tb3 
link to time dependent density functional theoriesc5L^. 
A overview about different approxirnation schemes are 
given inEl. We will see in paragraph III D that in order 
to fulfill the third order frequency sum rule we have to 
have an additional effective mass beyond £xc- 



required small wavevector limit of the so called screened 
structure fmiction (^9|) takes the form [11 ~ n"'j'^] 



> = hm ^ 



1 - e-^" 
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where fq{uj) = fq{Lo) — Vq. Now we observe that 
dp 
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vanishes for small q. Therefore we have to perform the 
limit in (MS) in the distribution sense to obtain 
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It is not difficult to see that the second part vanishes and 
we obtain 
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which agrees with ([l8|). 

Therefore, we have shown that the dynamical local 
field ( ^6|) from the response function (33) leads to the 
same compressibility (^9|). This gives besides the com- 
pressibility sum rule already checked a second proof that 
we have derived a dynamical local field which leads to the 
correct compressibility. The static limit will allow now to 
complete compressibility and third order frequency sum 
rule simultaneously. 



D. Frequency sum rules 



C. Compressibility sum rule 

From ( p6[ ) one sees immediately that the required 
form for the compressibility sum rule ( pO| ) appears, since 
n"'j'^(0) — no(0). Therefore we have derived a dynami- 
cal response function and a local field which shows in the 
static limit the correct compressibility (IT^). This com- 
pressibility formulae ( p^ ) can be checked alternatively by 
calculating explicitly the frequency integral in (|l9| ) . The 



Now that we have the response function (^3|) at hand 
we can proceed and proof the frequency sum rules ( [2^ ) 
and (|25| ) explicitly. First we expand the polarization 
functions for large frequencies. The superscript denotes 
which conservation laws are obeyed, density (n), energy 
(E) and current (j) respectively. We obtain 

n [uj) = 5- - I + 2E^r + - — ^ 7 
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We see that the current conservation repairs some defi- 
ances of the Mermin-Das polarization function, 11'^ which 
obeys only density conservation, in that the imaginary 
part shows a different frequency behavior 
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The last formulae corrects a misprint in formula (23) o 
This different behavior of the imaginary part is also re- 
flected in different expressions for the third order moment 
(o(l/ci;^)) or third order sum rule. 

From (p^) and (|4^)-(^0|) we read off the sum rules 
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which holds for each fl", 11"'^, n"'j^^ and x- con- 
trast to that we will see now that the third order sum 
rule gives different results for the inclusion of different 
conservations laws. Using the polarization function in- 
cluding densifry, energy and momentum conservation we 
obtain from (49) 



/ — c.3lmn"J'i=(c.) = 2E^ 

J TT TO^ 



nq 
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We remark that according to ([47|)-(|5(]|) the Mermin-Das 
polarization (A14) including only density conservation 



or even including additionally energy conservation (A13) 
would yield an additional — ng^/mr^ term which is an 
artefact. This is repaired by additionally taking into ac- 
count momentum conservation. 

Comparing ( p5| ) with ( p4[ ) we see that just the last 
terms are missing. In order to obtain this sum rule we 
have to use the response function ( ^3[ ) and not the po- 
larization function for which this sum rule is actually 
designed. With (|36|) and (|3l|) one gets 



X(a;)=n"'^^^(c.) + 



m^ 



1 



1 



(55) 



Consequently, the third order sum rule (|3: 
correctly if one sets 



is rendered 



no(0) 5n vno(O) 

= -Vq{l + i{q)) + ^E. 
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' A 



(56) 



By the requirement (|5g) we have a possibility to fulfill the 
third order sumrule exactly from the dynamical response 
as well as static local field model. 



E. Consequences on selfenergies 

Let us now work out what that means for our selfen- 
ergy parameterization A and m. From ( |56| ) we obtain a 
determining condition for the effective mass and energy 
shift [I{q) - 1 + /(<?)] 



V,I{q) 



1 /n2(0) q^\ dlnm 



2mn \Uo{0) 



9 In 7 



SA 

Sn 



AE. 



(57) 



Since we work with the effective mass and shift param- 
eterization of the quasiparticle energy the difference be- 
tween E and AE, vanishes but we keep it for com- 
pleteness further on. 

Applying the small wavevector limit ( C25| ) and (|3^) we 
see now from (p^) for Coulomb systems 



lim 



SA 3 1 51nm 



2n^Kn a In 7 



AE 



5 kj 



7 + o(g2). 
(58) 



Oppositely from (|57|), the large wavevector or small dis- 
tance limit ( C26 ) and (^6|) reads 



lim 

q — ^oo 



SA 
Sn 



1 E dlmn 
3 dlnn 



AE--{l-go)V, + o{l/q^) 
(59) 



where the last term on the right side vanishes for 
Coulomb potentials and persists only for potentials which 
range falls faster than Coulomb. 

If we assume homogeneous systems, where A and m 
becomes independent of the wavevector, the equations 
( p8| ) and ( p9| ) determines the quasiparticle shift as well 
as the effective mass via 



a In? 



5k' 



f 



dlnn 



SA 
Sn 



2 n'^Ko 



1E_ 

3 n2 



5fe 



^7 

/ 



9 1 
2 EKo 



1 



AE 



\il^go)Vg + oil/q^) (60) 



where the last term vanishes for Coulomb potentials. If 
we remember the relation between the selfenergy a and 
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the effective mass (|17p we can determine from (|6Cj) the 
thermal averaged selfenergy and wave function renormal- 
ization Z — (1 — d^a)^'^ . For the special case of zero 
temperature and neglecting the usually small dpO part, 
it reads 



A = ctI 



10^ V v'M) 

' 



, l{n') 



(61) 



with 7 from (||) such that the effective mass takes the 
form 



niQ exp < mq Vq / dn 



, l{n') 
p)[n') 



(62) 



The general expression for finite temperatures is given by 

Since the needed expressions / in (|2^) or 7 in (^) are 
functions of the structure function which is given itself 
again by the response function ( p^ ) we have the usual 
self-consistent procedure analogous to (||) first introduced 
by Singwi Sjolandera but with another G{q). Here we 
suggest to obtain the effective mass m and energy shift 
by (|7l) using /(g). 

Alternatively one might use the well known experimen- 
tal values of the static structure factor Sq and determine 
by this way the effective mass and energy shift. This def- 
inition of effective mass has the advantage that the third 
order and compressibility sum rule of the response func- 
tion will be rendered exactly. Therefore equation ( |60| ) is 
the second main result of this paper. 

Let us remark that if we would have no effective mass 
but a mere density dependent selfenergy a like in den- 
sity functional theories, the requirement of (|5|) and (|59| ) 
corresponding to the compressibility and third order fre- 
quency sum rule cannot be fulfilled SHH|tJ|taneously. This 
was remarked in detail in literaturePIFi. By including 
the effective mass we can resolve this puzzle here. 



IV. NUMERICAL RESULTS AND DISCUSSION 

As a test example we will now consider the unpolarized 
electron gas in three dimensions (3D) at zero tempera- 
ture. The density parameter is the usual Bruckner pa- 
rameter, defined as the ratio of interparticle distance to 
the Bohr radius, Vg = (3/47rn)^/'^/ao. First we will give 
a simple quasiparticle picture (QP), where the A and 
the effective mass are determined as density-dependent 
constants from MC data. In the second step, we will 
allow that the A and the effective mass depend on the 
wavevector. This will lead to the self-consistent quasi- 
particle picture (SQP). 



In order to calculate the quasiparticle parameter A 
and the effective mass, jwa-employ the results of Monte 
Carlo (MC) simulations ^EJ. For a parameterization of 
the MC data seeLj. Iru was discussed the difference of 
interacting and free kinetic energies, Ai? = -^{Eint— E). 
This difference is given as (5 = Eint / E — l\yy the MC data 
of Ceperley and Alder. In our quasiparticle picture one 
has Eint/E = mo/m, and the effective mass is given by 
m = mo/(l + S). This allows us to determine the needed 
derivative as 



d5 



dlnm 



d\nn '6{l + S)dr, 



(63) 



Ts 


70 
MC 


S 
MC 


MC 


1 - 
SQP 


5(0) 
QP 


PV4-AE 


1 


0.2567 


0.036 


0.7276 


0.7240 


0.6626 


0.6971 


2 




0.091 


(0.8627) 


0.8627 


0.7666 


0.8346 


3 


0.2722 




0.9078 


0.9201 


0.8456 




5 


0.2850 


0.292 


0.9768 


0.9627 


0.9351 


1.0558 


10 


0.3079 


0.619 


0.9976 


0.9733 


0.9613 


1.1231 



TABLE I. Tlie long wavelength limit of the local field fac- 
tor, 70, and the small raiiae value of the pair correlation ac- 
cording to the MC data o£j. The value in brackets is an inter- 
polation. The quasiparticle (QP), self-consistent quasiparticle 
(SQP) as well as Pathak-Vashishta value together with AE 
(PV+AE) is given for comparison. 

The difference AE is given by the large-wavelength 
limit of the local field 



\im G{q) 

q^O 



70 



kf 

= lim — 



1 



y„ \ 2mn 9 In n 



limJ-{VqI{q) 

dlnm SA 

Sn 



AE) 



^ ^ (64) 



which is presented by 7oE3, see table U. This allows us 



to determine the quasiparticle energy -j^, since the last 



line of ( |64D is just (|5^). The first model can be called 
improved Pathak/Vashishta scheme (PV-|-A£^) while the 
second one together with the effective mass (£^) estab- 
lishes the quasiparticle picture proposed here. We re- 
mark that, in order to realize a certain Bruckner param- 
eter , the quasiparticle picture must be calculated with 
rs{l + S), since all formulae work with effective mass lead- 
ing to rs/(l -t- S). 
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r =1 



r =10 



QP 

jMonte Carlo 

PV+AE 



q/k, 

FIG. 1. The static structure factor for 
from the (PV)-r4he (PV+AE)- p«d (QP)- 
to the MC dat<£3 as presented iufcl. 



= 1 and Ts = 10 
model compared 



With this determination of the QP parameters as well 
as (PV+Ai?) from MC data we can now compare the 
results of the iteration scheme fq ^ Sq ^ I{q) fq. 
For zero temperature all formulae simplifies and we give 
them here for convenience in dimensionless momentum 
K = q/kf, energy Q = w/ef and distance R = rkj/Ti. 
Then the dimensionless effective local field is given by 



fK = 



f2s^V'^ fI(K) 2 



4 
97r4 



86 



4(1 + (5) drs 



1 - 



no(o,x) 4 



-r (65) 



with 70 from (||) and 7 = -i / dK{SK - 1) from (|) 



and the spin degeneracy for electrons s = 2. The (PV) 
model would consists only in the term I{K) on the right 
hand side of (^) and the (PV+AE') model takes into ac- 
count the first line of (|65|). The (QP) picture finally takes 
all terms into account where we have used (|6|) or ( |57| ) 
and the requirement (|6j). The dimensionless function, 

T^SS) = ^fSm P^°"'='^ figure! is given by 



n2(o,x) _ 1 3 " ~ 47? (1 ~ I'^liqifl 




Provided we know the effective local field, Jk in (65) 
the static structure factor can be obtained from (noh as 



S 



K 



3 

47r 



d£t Im ■ 



i-(/A- + (|f;)i/3-^)no(f^,x) 



(67) 



where we have used the zer o te mperature dimensionless 
quantum polarization from (Bl) 



no(f7, K) = -l + A_ L. In . 



2K -Q + K^ 



4i^2 - (J] + 74-2)2 2K + n + K^ 
— In ■ 



8K^ 





-i- } f iorK <2and\n\ <\K^ -2K\ 

2 1 fo, „ 2K\ < \n\ < \K^ + 2K\ 



(68) 



With the help of the static structure factor we have the 
pair correlation function (O) 



00 

OR = 1 + ^ J dKK sin (KR) {Sk - 1) 



(69) 



fro m wh ich the required I{K) function reads according 
to (|C2i) 



00 

I{K)^-2 J ^ig^^l)j,(KR) 



(70) 



with the spherical Bessel function j2{x). This function 
now enters (|65|) closing the iteration. 




FIG. 2. The static local field for = 2 from the (PV)-r-the 
(PV+A_E)- and (QP)- model compared to the MC dataCJ. 

In table | we compare the small distance value of the 
pair correlation ( pT| ) of the (PV+AE') model with the 
(QP) model. We see that for more dense systems the 
(PV+AE) model leads to correlations which are too large 
while the (QP) model is lower than the MC values at 
higher r^. In figure |l| we compare the static structure 
factor of the two models with the MC data. We see that a 
difference occurs between the (PV+AE) and (QP) model 
at higher Ts- 
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S 1 

(5 



r.=5 



-- Gp„ 
• MC data 

- - Gpv+AE 



q/k, 

FIG. 3. The static local field for = 5 from the (PV)-p-jthe 
(PV+A_B)- and (QP)- model compared to the MC dataEj. 

The difference between the models becomes more ap- 
parent if we plot the local field factor as in figures gjij 
We see that the simple (PV) model underestimates the 
MC data though it satisfies the third order frequency sum 
rule. This result is improved by adding the Ai? read off 
from the large wavelength limit of MC data. Further 
improvement is achieved in the (QP) picture. At smaller 
densities, = 2, the local field is overestimated at higher 
wavevectors which leads to the deviation seen in table |. 




FIG. 4. The static local field for Vs = 10 from the (PV^ 
the (PV+A_B)- and (QP)- model compared to the MC data 



The (QP) curves show a small hump ai q — 2kf in 
contrast to the (PV-f AE') model. This comes from the 
which is plotted in figure ^. 



function w^tS 



cr 
o 



FIG. 5 




^ of ( p?! ) or explicitly for three di- 



no(0) 

mensions at zero temperature 



While in 3D this hump at q — 2ku is not much pro- 
nounced, it shows up in 2D systemsE^. With the help of 
the formulae in appendix ^ one can present also the 2D 
results which should be devoted to another paper. 

According to (^^ the functional form of 112/110 on 
the right hand side should be equal to I(k) which is a 
smooth function according to figures Therefore with 
the simple quasiparticle picture considered so far one can- 
not satisfy the third order sum rule for all q. In order to 
achieve this, we must allow A as well as the effective mass 
to carry a ^-dependence. Of course this leads to a self- 
consistent quasiparticle picture, since the energies under 
integration of the polarization function now change their 
dispersion. The iteration scheme is therefore enlarged 
to fq ^ Sq ^ I{q) A{q),m{q) fq according to 
(p7|). We call this iteration here self-consistent quasipar- 
ticle picture (SQP). It satisfies the third order frequency 
sum rule and the compressibility sum rule simultaneously 
for all q. The results describe the MC data in figures ||- 
^. Also, the small distance value of the pair correlation 
function is now in better agreement with the MC data as 
can be seen in table B. 



V. SUMMARY 

We have derived a response function in the quasiparti- 
cle picture where the correlations are parameterized by a 
density-dependent effective mass, energy shift and a re- 
laxation time respecting density, energy and momentum 
conservation. The dynamical response function can be 
given in the form of a modified RPA including a dynam- 
ical local field. This local field leads in the static limit 
to the correct compressibility. The effective mass and 
quasiparticle energy shift are proposed to be determined 
by the requirement of the third order sum rule. This al- 
lows in turn to satisfy the compressibility sum rule simul- 
taneously. Since the effective mass is now a function of 
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the structure factor, one might use experimentally known 
values or simulation results to construct a more realistic 
quasiparticle picture. The consequences on microscopic 
expressions of the selfenergy are presented. The thermal 
averaged selfenergy and wave function renormalization, 
i.e. the frequency derivative of the selfenergy, are linked 
to the pair correlation function at small distances. 

While all the expressions we derived are valid for finite 
temperatures, we have compared them as a test exam- 
ple with the Monte Carlo data for an electron gas at 
zero temperature. We find an improved description by 
the Pathak-Vashishta scheme accomplished by an energy 
shift derived from MC data. The best agreement with the 
data is achieved by constructing a nonlocal quasiparti- 
cle picture allowing for a wavevector-dependent effective 
mass and an energy shift in the form of a self-consistent 
quasiparticle picture. 
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APPENDIX A: DYNAMIC RESPONSE 
FUNCTION 



In the former paperBc^ it was derived formally the 
density, current and energy response {XiXJiXe} = 
X{1, 0, 0}"^ of an interacting quantum system 



6n \ 




q-SJ] = 


XJ 


SE / 


\XE 



yCKt ^ ^ 




(Al) 



to the external perturbation V'^^*' provided the density, 
momentum and energy are conserved. This has been 
achieved by linearizing the kinetic equation for the one- 
particle density operator p in relaxation time approxima- 
tion 



P' 



(A2) 



where the relaxation was considered with respect to the 
local density operator p^ °- or the corresponding local 
equilibrium distribution function 



Fip) = 



eoip-QjR, t))-pL{R,t) 
T{R,t) 



(A3) 



This local equilibrium is given by a local time dependent 
chemical potential /i, a local temperature T and a local 
mass motion momentum Q. These local quantities have 
been specified by the requirement that the expectation 
values for density, momentum and energy are the same as 



the expectation values performed with F which ensures 
conservation laws. 

The correlations are shared in the kinetic equation 
in such a way that the energy operator E parameter- 
izes the density dependent quasiparticle energy or varia- 
tion of the energy functional in the Landau liquid (mo- 
mentum dependent) or density functional sense (momen- 
tum independent) and the collision integral is approx- 
imated lat, a conserving relaxation time approximation. 
While idfj the general density, energy and momentum 
dependent form of such parameterization has been dis- 
cussed we want to consider now only a special case of an 
effective mass and rigid shift parameterization 



£ = -V( 



1 



2m(n 



■)V + A(n) 



(A4) 



such that the variation £(n) = £q + S£ reads [p = {pi 
P2)/2, q = pi -P2] 



5£ pi\5£\p2 >= {Vo + ViP^)Sniq) 



(A5) 



with 



Vo 



6n 

6_J_ 

Sn 2m 



V,- 



1 91nm 
2nm dhin 



(A6) 



The response matrix ( |Al| ) can be giwn in terms of the 
polarization matrix V, see Eq. (29) o£3, which is the re- 
sponse of the kinetic equation without the self-consistent 
quasiparticle energies S£. The response reads 



(A7) 



with the matrices simplified for our considered case (A5) 

(A8) 



/ 31^0 +^^^4 0^ 
V - .gpqVb + Vpq^4 



91 5pq 9e 
V 9^ ^?epq 9ee ■ 



(A9) 



and the correlation functions [e — p^/2m-|- A] are defined 
as 



30 (w) 



dp f(p+f)-F(p-f) 
(27r)3'^e(p+f)-e(p-f)-c^-i0" 



(AlO) 



P. We keep the matrix 
case in order to con- 



Explicit formulae are in appendix 
notation ofc3 also for this specia 
vince the reader about the technical usefulness of such 
notation. The 3x3 polarization matrix V = {n„„i} con- 
tains the corresponding density, momentum and energy 
polarizations as 



(All) 



Sn \ 






6q-j] 




^ind 


6E J 


Vnai/ 
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In fact as found in@ it is possible to express the density 
polarization function Hn = n"J'^ including density, cur- 
rent and energy conservation by a simpler one containing 
density and energy conservation H"'^ as 



1 



1 



1 



1 



n"J.E(^) nn^E(^) n"'j(tj) n"(tj) 



(A12) 



This shows that the momentum conservation leads sim- 
ply to a dynamical local field correction. The energy and 
density conserving polarization function reads explicitly 



n"-^(tj) = (1 - iut) 



gi(c^ + ^)gi(0) 
hi 



. (Mi(O)-M.(O))^ 
hi{h1-heehi) 



(A13) 



where we use the abbreviation /i^ = ^^(w -|- — 
LOTigfj,(Q). The first part is just the known Mer- 
min -Das pplaiization function including only density 
conservatiorOO 



no(o) 



(1 



ILOT 



gi(^ + ^)gi(0) 



(A14) 



Now we want to give the full density response function 
X — Xii according to ( Al). Due to the special considered 
case (AS) and consequently ( Ag), t he density response 
function x can be written from (A7) into [Xln = n"'j'^] 



x(w) 

where 

ni3(w) = 



UJ 



1 - VoTl'^'^'^iuj) - 2mF4ni3(w) 

n"'j''^ nftn2(o) - i™n2(w -t- ^) 

2m nhno(O) - iTUjUoiuj + 

n2(o)-no(o)n4(o) 



(A15) 



(A16) 



are e xpre ssed in terms of moments of the correlation func- 
tion (Bl). The response function ( |3^ ) is the main result 
of this paper since it gives the consistent response func- 
tion for the quasiparticle consisting of effective mass, en- 
ergy and relaxation time. 



APPENDIX B: EXPLICIT FORMULAE OF 
CORRELATION FUNCTIONS 



The different occurring correlation functions ( AlO) can 
be written in terms of moments of the usual Lindhard 
polarization function Ho 



dp FiP+l)~FiP-l) 



as 



9i = IIo 

El 

2m 

Ik 

2m 

pq = mujll2 



(Bl) 



5. 

9ep- 



5(pq)^ = -mq^n + m^tj^IIo. 



(B2) 



Here s is the spin degeneracy and D gives the dimension 
of the system. While all formulae in the text are written 
for the three dimensional case they hold equally for one 
- and two dimensions. 

For practical and numerical calculations we can rewrite 
the n„ by polynomial division into 



119 = —mn 



-Ho + n2 



-Bn 



q q 

2m^uj^~ ~ 14 2c.fl fori:> = 2,3 



q 



where the B^ are the projected moments perpendicular 
to q and read 



n2 - s 



dp 



■(p- 



pq ^2^(p + f)-^(p-f) 



/■ , ,„ r 2 for £> = 3 1 
= m ] d^Box|^ forZ? = l,2/ 



« toTBo X 

n4 = s 



2 forD = 3 



1 fori) = 1,2 

4 



dy , _ pq ^4^(p + §) - ^(p - §) 

{2tt)D g2 q^ 



m 

8 for Z) = 3 
3 forD = l,2 

8 for £1 = 3 
3 forD = l,2 



dfi' / dfi"UQ X 



(B4) 



The corresponding last identities are valid only for non- 
degenerate, Maxwellian, distributions with temperature 
T. The general form of the polarization functions is pre- 
sented as an integral over the chemical potential /i of the 
Lindhard polarization Bq. This is applicable also to the 
degenerate case. 
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1. Long wavelength expansion 

In real situations it is often helpftil to have the smaU 
wavevector expansion of the various occurring correlation 
functions. With the help of (B4) and (B3) this can be 
tremendously simplified if the expansion for XIq is written 



fr ( \ 32mg2p 



E2 

4 



64^ 



Uoiiu) 



bj J (27r) 
24m4w 



dp ( {q-pf ^, ^ q^{q-v? J,,, 



The static case (B6) yields with (B7) 

^2 



no(o) = -a^n 



' din 



12m 



(Bll) 



(B12) 



F'\+o{q') (B5) 



and with (gSj) and ( |B4| ) 



where i^' = etc. For the static case we have 



112(0) 



-2mn + —dfj_n 



no(o) 



{2Tr)D \ 8m 24m2 



16 



F" 



0{q'). 
(B6) 



n4(0) = -—m^E 

o 



2mq^ 



(B13) 



We give now the one, two and three dimensional case 
separately. 



a. 3D case 

Since dpF ~ —pdf^F/m partial integration gives 



and apphed to (B5) one gets 



E + oiq'^). 



(B8) 



During the text we use also the small and large wavevec- 
tor limit of the static 3-D polarization functions. Since 



no(o) 



-An^-^E + o{l/q^) 



we obtain from (B3) and (B4) 



n2(o) , 

' —mn 
which leads to 

n2(0) 
no(o) 



— 3mn + o{q^) 

16 TP 12Sm^ rp i „/i /„6 



oil/q') 



2Em 
3n 



(B14) 



(B15) 



(B16) 



The density, n, and energy, E, and higher moments read 
in terms of 



as 



1 



T{n) J e'^-^t^ + 1 




^,3 = 



2to^ 

2 \ 3 



15 s 



4/32 A3 ■ 



(B9) 



(BIO) 



6. 2D case 



For two dimensions we have instead of (B7) 



dp 
(2^ 



■F'G{p) = TO 



dp 



■FdpGip) 



which applied to (B5) yields 



no(w) 



27rg2 



Girq 



■E 



and from (B3) and (B4) 



n2(w) - ^E+-^E2 



mu!^ 
,2^ 1C^„4 



U^iuj) = —^E2 + ^E3. 

The different occurring moments reads here 



(B17) 



(B18) 



(B19) 



With the help of (B4) and (B3) one writes down imme- 
diately the higher order correlation functions as 
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n= (1) 



E = 



^2rn 



^2 = 



— 



P 

2m 
2to 



(B20) 



The static case is now analogously and reads with (B17) 
from (^) 



no(0) = -a,n+^9> 
2 



112(0) = -mn+ I^^m" 



04(0) = -3m^E 



mq 



(B21) 



c. i-D case 



For one dimensions we have instead of (B7) 



'^P F'G{p) ^mf ^Fd,{^: 



{2n 



(B22) 



which applied to ( |B5D yields 



? 4 
llo(cLij = -n — — — - — -OuTi 



6g4 



?1 + 7 jE2 



E 



2a^ 

U2{lj) = -^E „,„ , . 

n4(^) = ^r^S2 - + 7tA^3. 



3u;2 

The static case is 

no(o) = -d^ 



25^4 



(B23) 



o,,n 



mq 



12m 

112(0) = -mn + j^df,n 
114(0) = - = -ern^i; 



The occurring moments reads here 
2 



(B24) 



9 M 

, 3 



d^l'E 



3 s 



I6/32 A 



h/2 



E:, 



P 

2m 



dfi'E2 



15 s 



32/33 ;v 



/7/2- 



(B25) 



2. Large frequency limit 

The large frequency limit can be given analogously to 
the forgoing chapter. We restrict here to give the expan- 
sion for 3D 



2q* 



-E^ 



q 



Am^uj^ 



n + o{—). (B26) 



Please re mar k the difference to the large wavelength ex- 
pansion (B8). The corresponding higher order correla- 
tion functions are completely analogousl y giv en by the 
methods of the foregoing chapter. With (B4) and (B3) 
one gets 



n2(w) 
n4(w) 



E- 



Aq^ 

32to9^ 
15^2 



E2 



1 



Zm^oj' 



-E + o{—A 



Eo 



64g4 
35^4 



E. 



8q^ 



1 



E2+0{-r). 



(B27) 



APPENDIX C: PERTURBATION THEORY AND 
FREQUENCY SUM RULES FOR 1,2,3 
DIMENSIONS 

The external potential is adiabatically switched on 

y<='^t(r,t) = y(r)e°*e(-t) (CI) 

and induces a time dependent change in the Hamilton 
operator 



SHit) 



drn{r, t)V''''\r, t). 



(C2) 



The variation of the density matrix operator p{t) — 
p-\-5p{t) can be found from the linearized van-Neumann 
equation as 



Spit) 



[SH,po] 



(C3) 



where it has been assumed that the perturbation is 
conserving symmetries of the equilibrium Hamiltonian 
[Ho,6p]=0.^ 

The variation of the density expectation value 5n = 
TrSpn is consequently 
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Sn{r,t)=i / dt' dr'V{r',t'){[n{r,t),n{r',t')]). 



(C4) 

Since in equilibrium the commutator is only dependent 
on the difference of coordinates and times we can define 

x{[n{r,t),n{r\t')]) (C5) 



from which we obtain the Fourier transform of (C4) to 
Sni,,.) = V-\,,.) f'^^^^^ (C6) 

J TT U ^ UJ — IV 

where V'^^^{q,Lj) — V°^^{q)/{0 + iuj). This is of course 
identical with (nl). 



1. Sum rules 

Now one can derive the first and second order sum 
rules of the imaginary part of the response function (|C5| ). 
Therefore we generalize the definition ( p5| ) to nonequi- 
librium and finite systems 

Inix((7,t^,i?,i) = j dre''^^ j dre"*'?'' 

X {[n{R + r/2, t + r/2), n{R - r/2, t ~ r/2)]). (C7) 

Higher order mcaiients can be expressed by correlation 
functions as welH. Here we restrict to the lowest two 
orders and rederive it in conventional way. 

If we assume further equilibrium but finite systems we 
can define an avera ged response by applying spatial av- 
eraging / dR/V to (C7) such that we obtain 



Imx(<?,^)=-^ j dTe^^'{[n{q,t)M~qM)- (C8) 

From this expression it is easy to see that the first two 
frequency sum rules read 

— wlmx(g,t^) = ~^{[idth{q,t)\t^a,n{~q,Gi)]) 

TT V 

—u;Hmx{q,Uj) =-h[{^^t)^n{q, t)\^, h{-q, 0)]). 

TT V 

(C9) 

Using the Heisenberg equation idtft = [n, H] and 
H 



dp 



(27r)^ 2m''f 



1 f dpdpidp2^. . , . . , 



2 J {2n)^^ ''p"pi"p2"P2+p"pi-p 
_dp_. + . 



(CIO) 



we can express the sum rules (Oh as 



duj 



1 



— ujlmx{q,oj) = 77 (ki 



TT 

duj 



V 



1 



{[qj,.,[qj^„H]]) (Cll) 



TT ■ V 

where the divergence of the current operator reads 

q 



qjq 



dp 2pq 



-a,p Clp+q. 



{2tt)D 2m 

Performing the last commutators one obtains finally 



(C12) 



dhJ 



ujlmx{q,uj) 



and 



J^u;^lmx{q,u;) 



nq" 



{nq=o) 
m V 



Sq^ f dp 



(C13) 



b-'?)^ (ftp Op) 



1 



4to3 mWj{2TT)D 

f_^P_T. // ^2 fdpidp2 . 
mW 1 {2tt)d Py^'P''^' y (■27r)20°'PiS2"P2+p+?api-p 

\2D '^Pi '^P2^P2-p'^Pi+'P 



-{{p-qf~q\p-q)) 



{2nY 



(C14) 



Since we had symmetric expressions 4>(j>) = (f'i^p) the 
second term leads just to kinetic energy density 



dp 



31 

m?Y J {2tt)D 

2q'^ f dp p^ 



mP-V J {2tt)D 2m 



{p-qf{dpdj,) 

{dp dp) X 



3 for ID 
1 for2,3Z) 



2q^E 



3 for ID 
1 for 2,31) 



(C15) 



Now we are going to express the last 4 creation and 
annihilation operators by the structure function itself. 
Therefore we use the definition of the pair correlation 
function 



dj^dj^dr^dr^) = grn{ri)n{r2). 



(C16) 



Applying the spatial averaging j dR/V where R = {ri + 
r2)/2 and Fourier-transform the difference ri — r2 into q 
we obtain 

f dpidp2 . 

= J dre-'^-'^grJ ^n{R + r /2)n{R - r/2) 

= n^ J dre-'''"{gr - 1) + n^{2TT f5{q) 

= n{Sci~l) + ji^{2TTf5{q) (C17) 

where we neglected spatial gradients in the density and 
used (|ll|) for the last step. Using ( |C17| ) in ( |C14| ) we 
obtain finally 
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-Lu^l-mx{q,uj) = 



nq" 
Avrfi 



2q^E 



n 
m 



^q%{liq) + 



3 for ID 
1 for 2, 31? 

ii{q) for ID 
for 2, 3D 

(C18) 



with (M) 



dk 



n J (27r) 



(C19) 



where we understand Sk — (27r)'^5(fc). For the one- 
dimensional case an extra term appears 



h{q) = - / jTr-ASk-l + nSkfl)-^. 
n J [2tt) qVq 



(C20) 



In the following we restrict to the required formulae for 
the three dimensional case. One can Fourier transform 

i{q) = - f drgr{l - cos {q ■ r))^^^^ (C21) 

J q ''q 



which was first given by Pufia. This correct form lead s 
unavoidably to the appearance of the Sp terms in ( CIS ) 
very often overseen in later papers. For asymptotic ex- 
pansions, however, we have to be careful that — 1 is the 
object which re nders spatial integrals finite. Therefore 
the Sp terms in (C19) have to be considered separately 



1 



dk 



n J (27r)' 



(fc • qf Vk 
q* Vq 



= -1 



(C22) 



such that we obtain instead of (C21) 



dr{gr - 1)(1 - cos {q ■ r)) 



{q ■ drfVr 
q^Vg 



(C23) 



For Coulomb potentials we can further simplify 

oo 

dr 



Hi) 



-{gr " i)j2{qr) - 1 



(C24) 



with the spherical Bessel function j2{x). From this ex- 
pression one sees the small wavevector limit 



/(9) = -^ I drr{Sr-l)-l + o(q') 



||J.-l + o(,') 



(C25) 



where we have used ( pi] ) and the definition . The long 
wavevector limit takes the form 



f dx 

I{q) = -2 J —{g.iq " l)i2(x) - 1 



= -2(30-1) 1 ^j2(x)-l + o(l/g2) 



= ^(l-5o)- 1 + 0(1/?'). (C26) 
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